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We examine the absorption spectrum of electromagnetic radiation from excitons, where an exciton 
in the Is state absorbs a photon and makes a transition to the 2p state. We demonstrate that 
the absorption spectrum depends strongly on the quantum degeneracy of the exciton gas, and that 
it will generally manifest many-body effects. Based on our results we propose that absorption of 
infrared radiation could resolve recent contradictory experimental results on excitons in CU2O. 
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The phenomenon of Bose-Einstein condensation has 
attracted much attention in the recent years, and many 
experimental groups have reported the formation of Bose- 
Einstein condensates in vapors of alkali- metal atoms [jj] . 
Excitons, bound states of electrons and holes in semicon- 
ducting materials, are another candidate for undergoing 
this phase transition pj. Since excitons are composite 
particles consisting of two fermions, they are expected 
to obey Bose-Einstein statistics in the limit where their 
thermal de Broglie wavelength becomes comparable to 
their interparticle spacing, provided that this spacing is 
much larger than the exciton Bohr radius. 

Excitons in CU2O, which is a dipole forbidden mate- 
rial, have some important advantages in this respect || . 
Their radiative lifetime is long and their binding energy 
is much larger than the typical thermal energies. In addi- 
tion the effective electron and hole masses are isotropic, 
and finally there are no bound states between excitons, 
i.e., biexcitons, or an electron-hole-liquid phase. 

A lot of effort has been made in order to create a Bose- 
Einstein condensate of excitons in CU2O 0||. To deter- 
mine the density and the temperature in the above ex- 
periments, the phonon-assisted recombination spectrum 
was fitted to a Bose-Einstein distribution, which gave the 
chemical potential and the temperature - two essentially 
independent parameters. Given the total exciton mass, 
the density was then evaluated to be on the order of 10 18 
cm~ 3 , while the temperature was on the order of 20 - 30 
K, higher than the lattice temperature which was kept 
at about 5 K. Therefore, according to this approach, the 
exciton gas was very close to the phase boundary for 
Bose-Einstein condensation, and the angular-momentum 
singlet-state (para)excitons were the species reported to 
actually cross the phase boundary. 

More recent experiments || have, however, questioned 
the older method of determining the density and the tem- 
perature. In these experiments the number of excitons 
was determined directly and, with a relatively reliable es- 
timate of the volume of the exciton cloud, the density was 
found to be two to three orders of magnitude lower, i.e., 
around 10 16 cm~ 3 , where the exciton gas should show no 
sign of quantum degeneracy. 

One, therefore, needs to find a reliable method of deter- 
mining the density and in particular the degree of quan- 



tum degeneracy of excitons. In this study we propose 
that measuring the absorption spectrum of infrared radi- 
ation which induces transitions of the excitons from the 
Is to the 2p state can indeed resolve the discrepancy. 
Our study is also directly applicable to other systems, 
like excitons in quantum wells, and thus it could help 
resolve some other experimental observations which are 
controversial 0. 

Oktel and Levitov fi have examined a similar prob- 
lem as the one we consider here, in the context of op- 
tical excitations of hydrogen atoms and have studied 
the many-body effects that show up in the absorption 
spectrum, for an effective contact potential between the 
atoms. Our approach is equivalent to theirs in the limit 
of equal masses between the excitons in the Is and the 
2p state. In another study Pethick and Stoof || have 
considered a more general form of the interatomic poten- 
tial. Finally Lewenstein and You (Toj| have examined the 
enhancement of scattering of light from a gas of bosonic 
atoms as they form a Bose-Einstein condensate, and have 
suggested that this effect could be used for the detection 
of the condensate. 

This Letter is organized as follows: We first consider 
the case of an ideal exciton gas, and examine the rele- 
vant energy scales that enter the problem and also derive 
simple expressions for the absorption spectrum. We then 
examine the problem of an interacting exciton gas within 
the Hartree-Fock approximation and find that the inter- 
actions can have a very drastic effect on the absorption 
spectrum. Finally we present our results with the inter- 
actions included, and show that infrared absorption can 
be used in order to determine the degree of quantum de- 
generacy of excitons, thus proposing an experiment which 
could resolve this issue. 

Consider the process in which an exciton in the Is 
state, with momentum Tik, absorbs a photon of momen- 
tum hc{, making a transition to the 2p state with momen- 
tum ft(k + q). The conservation of energy in this process 
implies that 

ei' + ftcg = e& q , (1) 

where = Ei + h 2 k 2 /2mi with Ei being the binding 
energy of the i state, and rrii being the total exciton 
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mass in the state i. In C112O mi, ~ 3m, where m is 
the electron mass, is larger than the sum of the effective 
electron and hole masses as a result of the small Bohr 
radius of the Is excitons, a \? ~ 5.3 A compared to the 
lattice constant a; « 4.26 A On the other hand, the 
Bohr radius of excitons in the 2p state a 2 g is given by 
the hydrogenic formula which yields ~ 44 A jO]. Since 
a 2 g ^ ai, m,2 P is expected to be equal to the sum of 
the effective electron and hole masses, which is w 1.68m. 
In Eq. (Q) there are two distinct energy scales, i.e., the 
energy separation AE = E 2p — E\ s ss 55 meV, and the 
thermal energy h 2 k 2 /2m ~ ksT, which is of order 1-10 
meV. Since Tiqc ~ AE and h 2 k 2 /2m ~ fcsT, we get that 
q/k ~ AE/y/mc 2 k B T ~ 1CT 3 . Therefore h 2 qk/2m ~ 
10" 3 /fc B r, and h 2 q 2 /2m ~ 10" 6 fc B T, which allows us to 



neglect the corresponding terms in Eq. (|lj) 
terms of k 2 (q), we obtain 
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Solving in 
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which gives the magnitude of the momentum Kk of 
the exciton in the Is state that absorbs a photon with 
wavevector q and is excited to the 2p state. In this 
approximation, for hcq = AE only the excitons with 
k = can participate in the process; however for a Bosc- 
Einstein condensate there is a macroscopic number of 
excitons with k = 0, and therefore the absorption spec- 
trum has a pronounced peak, with a strong temperature 
dependence. To see this more clearly, let us calculate 
the rate of this dipole-allowed process of absorption of a 
photon. With the approximate expression for the conser- 
vation of energy of Eq. (||) , the rate Tt of non-condensed 
excitons in the Is state absorbing a photon and making 
the transition to the 2p state is given by 
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M q |X s (l + ™ k P +q )/ q 6(hqc - AE*), (3) 



k 



where M q is the matrix element of this process, is the 
distribution function of species i (Is or 2p excitons), / q 
is the distribution function of the incoming photons, and 
Ai?k = e k p — e lt ■ Neglecting the occupation number of 
excitons in the 2p state, «k+q ^ 1> f° r monochromatic 
radiation Eq. (ft) implies that 



T T oc (hcq - AE) 1/2 nl s o d(hcq - AE), 



(4) 



where 0(x) is the Heaviside step function, and the magni- 
tude of ko is given by Eq. (||) . The above result expresses 
the fact that the absorption spectrum is proportional to 
the density of states times the distribution function calcu- 
lated at a wavevector with a magnitude given by Eq. @ . 

For a Bose-Einstein condensed exciton gas with Nc 
excitons occupying the k = state, the rate Tc of the 
same process is simply 

T C = ^-iM^Ncil + n^f^Sihqc-AE), (5) 



or Tc oc Nc 8(hcq— AE). Therefore the absorption spec- 
trum (which is proportional to the decay rate) of an ideal 
Bose-Einstein condensed gas has a strong peak with a 
height that scales as Nc- However, as we show below, 
the interactions can modify this picture drastically. 

We thus turn to the more realistic problem of an in- 
teracting Bose gas. We start with the Hamiltonian H 

H 

H = 1^ e k < a k + L a k+q a k'-q a k'a k 



k,k',q 



£4 p ^ k + ^ E *i + A-^ (6) 



V 



k,k',q 



where V is the volume of the gas. In the above Hamilto- 
nian we assume that the excitons interact with an effec- 
tive contact potential, with Uij = 2nh 2 aij / jiij being the 
strength of the effective two-body interaction. Here ay 
is the scattering length for collisions between excitons in 
the states i and j (1 for the ls state, and 2 for the 2p 
state). The reduced mass ^ entering the above expres- 
sion is given by /Xy = m,mj / '(m,; + mj). Also a k (6 k ) 
and 0,^(6^) are annihilation and creation operators for 
an exciton with momentum k in the ls(2p) state. 

Let us now consider the ground state of the sys- 
tem with N excitons, which we denote as |0) = 
|iV"c, AT kl , AT k2 , . . . , Nk e ,...), where iV kj is the occupancy 
of a state with momentum k^ . Initially we take all the ex- 
citons to be in the ls state. Since we consider a Bose gas 
both in the normal, as well as in the condensed regime, we 
assume that there is one state that can get populated by 
a macroscopic number of excitons, Nc, and thus Nc can 
get as high as the total number of excitons, N, whereas 
the iV ki are of order unity. 

We now examine such a system when one creates ex- 
citations of the excitons from the ls to the 2p state with 
the action of some laser pulse. If an exciton with mo- 
mentum k e is excited to the 2p state with momentum 
k e = k e + q, where q is the wavevector of the laser 
light, since q is very small, we shall assume that we have 
vertical transitions, i.e., q = 0. We denote the excited 
states as |$ oxc ,k e ) = |k e ; Nc,N kl , N kz , ■ ■ -,N ke - 1,.. .), 
which are the basis vectors of our problem. The num- 
ber of such states is N — Nc + 1 = Nt + 1, where 
Nt = N — Nc is the number of ls excitons in states 
with k ^ 0. The laser beam that excites the excitons 
from the ls to the 2p state creates a superposition of the 
states I'&excke) jjJI- Thus, in order to determine the ab- 
sorption spectrum, we consider the matrix with elements 
H itj = (SexckJ-Hl^exck,) - (0|J3"|0>. One finds that 

H id = S i>:j [AE ki + f7u(n ki - 2n) + U 12 (n - n k J] 

+ Uu^n^n^, (7) 

where n = N/V and n ki = n k s = N-^JV is the Bose- 
Einstein distribution for the ls excitons. Let be the 
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components of an eigenvector with eigenvalue E. Start- 

and sum over 



ing from the eigenvalue equation ^ 
we solve in terms of multiply by 
i. Eliminating the factor J^j 
ing equation, the eigenvalues of Hi 3 
the roots of g(E) — 1 = 0, where 




the paraexcitons is ~ 10 s. The radiative lifetime 
in the 2p state rf p is smaller by a factor of (fc 7 a^) 2 , 
since the transition is dipole allowed, where fc 7 is the 
from the result- wavevector of the emitted photon. Since fc 7 = E g /hc, 
are then given by where E g w 2.17 eV is the gap energy, fc 7 sa 10" 3 A" 1 . 



Therefore t*p ~ 10" 



8 s, while Tp P 



10 _b s. The exciton- 



9(E) = E 

i=0 



Ui 2 n^ 



(8) 

Distinguishing the condensed state (i — 0) from the other 
states (i ^ 0), Eq. (j|) takes the following form in the 
thermodynamic limit 



9(E) 



phonon scattering time is on the order of 10 s [|13|. Fi- 
nally for a density as high as 10 18 cm" 3 , the exciton- 
E - [A£; ki + t/u(n ki - 2n) + U 12 (n - n kl )] ' exciton scattering time is on the order of 10" n s @, 

which turns out to be the shortest possible timescale, 
giving an energy broadening of less than 0.1 meV. 

We now analyze the results shown in Fig. 1. To pro- 
duce these graphs, we made use of the results of Ref. 
|]l5| , that an = 2.1a^ for paraexcitons, and assumed 
that an = 10 A. For the value of a\ 2 very little is known 
and for this reason we have considered both the case of 
positive (left column in Fig. 1), as well as negative (right 
column in Fig. 1) values. The ratio |ai2/an| is expected 
to be larger than 1, since a 2 ^ = 44 A, which is much 
larger than « 5 A. We made the conservative choice 
|ai2/an | = 2, although this ratio could be larger. We 
also considered a temperature of 10 K for the exciton gas 
in all the cases, and we varied the density from 10 16 cm" 3 
to 5 x 10 18 cm" 3 . With these values ksT ~ 1 meV, while 



U 12 n c 



E-[AE + U n (n c - 2n) + U 12 (n - n c )] 
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E - [A£ ki + n(U 12 - 2Un) 
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where nc — Nc/N. In the above equation there are in 
general three limiting cases, depending on the ratio of 
the interaction energy nU\ 2 , to the typical kinetic en- 
ergy A-E ki , which is on the order of the thermal energy, 
ksT. In the limit nU± 2 <C ksT, one recovers the results 
we found earlier for the ideal Bose gas. In the oppo- 
site limit, nU\ 2 3> fcgT, the behaviour of the system of 
excitons is "collective" . A graphical solution of the eigen- 
value equation shows that in the condensed phase, where 
both Nc and Nt are of order N, there are two strong 
modes, which give rise to two peaks in the absorption 
spectrum. There are also Nt — 1 solutions, which form 
a continuum corresponding to single-particle excitations 
of the thermal excitons. In the same limit nUi 2 3> ksT 
for a fully condensed gas as well as for a gas in the nor- 
mal state, there is only one mode, since then one has a 
one-component system. Finally when nUi 2 ~ fcgT the 
system behaves in a "mixed" way. In addition, the limit 
|wi s /m2p — 1| <C 1, is equivalent to the case nU\ 2 ^> ksT 
and Eq. (^) reduces to a quadratic algebraic equation, 
which gives the same result as the one derived by Oktel 
and Levitov in Ref. ||. 

By adding a small imaginary part in g{E) 1 i.e., g(E + 
in), where finite n results in homogeneous broadening of 
the energy levels, we calculate the corresponding imagi- 
nary part of the susceptibility [g(E + in) — l]" 1 obtaining 
the absorption spectra shown in Fig. 1. Broadening can 
be calculated from first principles |9j - however in the 
present study we assume small homogeneous broadening, 
choosing n — 10" 2 meV to produce the graphs in Fig. 1. 
The broadening of the energy levels is expected to be 
small, and this can be seen by examining the three basic 
mechanisms which contribute to that, i.e., the exciton- 
exciton elastic collisions, their scattering with the lattice, 
and their radiative lifetime. The radiative lifetime of the 
orthoexcitons in the Is state rl s is ~ 10 _5 s 101 and that of 



\nUi 2 \ is ~ 10 meV for n = 10 b cm" 3 , and ~ 5 meV 
for n = 5 x 10 18 cm" 3 . Figures l(a+) and (a-) show 
a completely classical gas, and since nU\ 2 <C fcsT, the 
system behaves like an ideal gas. In Figs. l(b+) and (b-) 
the excitons are essentially on the phase boundary for 
condensation and since nU\ 2 ~ ksE ', the system is in the 
"mixed" state where both collective and single-particle- 
like behaviours show up. Figure l(b-) shows these two 
distinct types of excitation, while Fig. l(b+) does not, be- 
cause the collective mode is buried inside the continuum. 
In Figs. l(c+) and (c-) the excitons are in the condensed 
phase with Nc/N w 0.48. This is the source of the sud- 
den appearance of the peak in Fig. l(c+). In Fig. l(c-) in 
addition to the two peaks, there is a contribution from the 
continuum that is hardly visible. Finally in Figs. l(d+) 
and (d-) Nc/N » 0.79, and since nU\ 2 is about 5fcsT, 
the spectrum is dominated by the collective behaviour, 
as the two peaks indicate. However, we remark that the 
Hartree-Fock approximation does not capture effects due 
to condensate fluctuations which may be relevant in the 
regime n\J\ 2 3> ksT. 

Let us now examine the possible experiment which 
could be performed in order for these effects to be ex- 
plored. The energy of the absorbed radiation would have 
to be in the infrared, with an energy of order AE w 55 
meV. The corresponding wavelength is about 20 /im, and 
it is comparable to the size of the cloud. Free-electron 
lasers provide tunable radiation in this regime. It is im- 
portant to mention that at such low energies the crystal 
is transparent and the absorption of radiation due to the 
process we study should be the dominant mechanism. 
Our analysis requires that the infrared pulse should be 
sufficiently long, so that its energy spread is much less 



3 




1 2 3 
Energy (meV) 





10 18 cm" 3 




(b+) 


J 





-8 -6 -4 -2 


2 


Energy (meV) 




10 18 cm" 3 




(b-) 


. J 






-4 -2 
Energy (meV) 





2x1 18 cm" 3 


J 





-4 -2 
Energy (meV) 



5x10 cm 
(*■) 



2 3 4 
Energy (meV) 



-4 -2 
Energy (meV) 



FIG. 1. The absorption spectrum on a linear scale, as a 
function of the energy of the photon that is absorbed. The 
absorption is in arbitrary units, but the scale is the same in all 
the figures. The energy is measured in units of meV, and its 
zero is measured with respect to AE. The scattering length 
an = 10 A, the exciton temperature is 10 K, and the density 
is 10 16 , 10 18 , 2 x 10 18 , and 5 x 10 18 cm" 3 from top to bottom. 
For the graphs on the left a\2 = 20 A, and for the ones on 
the right, a\2 = —20 A. If fj, is the chemical potential of the 
gas, -fi/fef = 3.7 in (a±), and 5 x 10" 4 in (b±), while 
N c /N = 0.48 in (c±), and 0.79 in (d±). 



than the energy width of the structures shown in Fig. 1. 
An advantage of the method we suggest is that it pro- 
vides an independent method of probing the kinetic en- 
ergy distribution of excitons. The difference between the 
uppermost graphs in Fig. 1 and the lowest is pronounced, 
and one should be able to distinguish clearly the degree 
of degeneracy of the excitons. In addition, this method 
does not depend on the strength of the phonon-assisted 
recombination line of paraexcitons, which is very weak, 
and since it is close to other much stronger lines, observ- 
ing this line is very hard §-§. 



In conclusion, we have demonstrated that the absorp- 
tion spectrum of electromagnetic radiation which induces 
transitions of the excitons from the Is to the 2p state is 
strongly affected by the quantum degeneracy of the ex- 
citon gas and by many-body effects. We have thus con- 
cluded that the absorption spectrum can provide a pow- 
erful tool for resolving recent experimental contradictory 
results on excitons in CU2O, and in quantum wells. 
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